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ABSTRACT 

Accretion onto compact objects plays a central role in high energy astrophysics. In 
these environments, both general relativistic and plasma effects may have significant 
impacts upon the propagation of photons. We present a full general relativistic magne- 
toionic theory, capable of tracing rays in the geometric optics approximation through 
a magnetised plasma in the vicinity of a compact object. We consider both the cold 
and warm, ion and pair plasmas. When plasma effects become large the two plasma 
eignemodes follow different ray trajectories resulting in a large observable polarisation. 
This has implications for accreting systems ranging from pulsars and X-ray binaries 
to AGN. 
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1 INTRODUCTION 

A considerable amount of effort has been invested in at- 
tempting to reproduce the spectral properties of accret- 
ing compact objects. A great deal of this work has been 
concerned with fitting the unpolarised flux with an un- 
derlying, phv^K2aUv_jnotivatedrnoo^ of an accretion flow 
(see e.g. Blandford fc Beeelmanlll999l : iN^rxav an fc Yi|ll994t 
iQuataert fe Gruzinovl 120001) . These models have met with 
some success, even being able to make testa ble predictions 
regar ding the accretion environment (see e.g. lNaravan et all 
1998). Because these models are primarily concerned with 
the physical structure of the accretion flow, they ignore the 
effects that the combination of dispersion and general rela- 
tivity will have upon the spectra. Far from the compact ob- 
ject this may not matter (r S> M). However, for the emission 
originating from near the compact object, this combination 
can be crucial. 

General relativistic vacuum propagation effects have 
been extensively studied in both the polarised and unpo- 
larised cases. In some systems gravitational lensing has been 
shown to have detectable c onsequences in cert ain regions of 
the spectrum. For example. iFalcke et al 1 fcooofr have argued 
that the black hole in the Galactic centre may be imaged 
directly at millimeter wavelengths as a result of gravita- 
tional lensing. In addition, general relativity has been shown 
to have a depolarising influen ce upon photons passing near 
the compact object (see e.g. ILaor et all ll99Ct lAgolll99*7i 
IConnors et"ai]ll980l) . However, the fact that these studies 
ignore plasma effects make them inapplicable to thick disks 
and at frequencies near the plasma and/or cyclotron fre- 
quencies. 

On the other hand, astrophysical plasma effects have 
also been studied in detail, although primarily in the con- 
text of nondispersive propagation effects upon the po- 



larisation, e.g. Faraday ro t ation and conversion (see e.g. 
Sazonov fc Tsvtovichl Il96£t ISazonovl ll96St IJones fc O'Dell 
1977allbt iRuszkowski fc Beeelmanl |2002|) . Weak dispersion 
has been considered in the form of scintillation (see 
e.g. iMacauart fc Melrosel l2000l) . While this can lead to a 
high degree of polarisation variability, it has a vanishing 
time/spatially averaged value and doesn't otherwise affect 
spectral properties. In contrast, strongly dispersive plasma 
effects have been extensively studied in the context of radio 
wave s propagating through the ionosphere (see e.g. iBuddenl 
Il964l) . Here it has been found that dispersive plasma effects 
can play an important role in determining the intensity and 
limiting polarisation of the radio waves. None the less, nei- 
ther of these type of plasma effects have been studied in 
conjunction with general relativistic effects. 

There have been some attempts at treating both gen- 
eral relativistic and plasma effects. However, these have 
been restricted to eit her unmagnetised plasmas (see e.g. 
iKulsrud fc Loeblll992l) or to nondispersive emission effects 
(see e.g. lBTomlev^'t^l]l200 j) . Both of these have only lim- 
ited applicability for realistic accretion flows. 

We present a fully general relativistic magnetoionic the- 
ory. This is a natural extension of the previous work combin- 
ing both general relativistic and plasma effects upon wave 
propagation in the geometrical optics limit. This will be pre- 
sented in four sections with [JH developing the theory, Sj3] 
presenting some simple example applications, and con ~ 

taining conclusions. Throughout this paper the ( h H — h) 

metric signature will be used, and h = c = 1. 

In a subsequent paper II we will discuss the details of 
performing radiative transfer in general relativistic plasma 
environments. These have been expressly neglected here in 
the interest of clarity. 
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2 THEORY 

The natural place to begin a study of plasma modes is 
the covariant form ulation of Maxwell's equations (see e.g. 
iMisner et aT1ll973l) : 

V^F"" = 4tt J" and V/F" M = , (1) 

where F v>1 = VA* 1 - V M A" is the electromagnetic field ten- 
sor, *F"'< = \z"^F a p is the dual to F M " (e^ a/3 is the Levi- 
Civita pseudo tensor) , and J v is the current four-vector. In 
order to close this set of equations, a relation between the 
current and the electromagnetic fields is required. For the 
field strengths of interest here, this will take the form of 
Ohm's Law: 

r = a\F^u a , (2) 

where uf 1 is the average plasma four-velocity and a"^ is the 
covariant generalisation of the conductivity tensor, defined 
by this relationship. As a result of the anti-symmetry of F^ v ', 
the conductivity will in general have only nine physically 
meaningful components, namely the spatial components in 
the slicing orthogonal to w M . Nonetheless, in order to investi- 
gate the behaviours of plasma modes in a general relativistic 
environment, it is necessary to express the conductivity in 
this covariant fashion. 

This can be more naturally expressed in terms of = 
F^Uv and = *F M "it„, the four- vectors coincident with 
the electric and magnetic field vectors in the locally flat 
centre-of-mass rest (LFCR) frame of the plasma. In terms 
of E^ and B^, the electromagnetic field tensor and its dual 
take the forms 

F" v = u»E v - E^u" + e^ a& u a B l3 , (3) 
= u"B v - B^u v + u a Ep . (4) 

Inserting these and Ohm's law into Maxwell's equations 
yields eight partial differential equations, 

V M (ifE* - E v u? + e"» a& u a Bp) = ^a^E" , (5) 

V M (u"B» - B v lf + e v ^ u a E^=0, (6) 

which may be solved for E^ and _B M given an explicit form 
of the conductivity. 

2.1 Geometric Optics Approximation 

The general case can be prohibitively difficult to solve for 
physically interesting plasmas. Fortunately, the problem can 
be significantly simplified by making use of a two length 
scale expansion (also known as the WKB, Eikonal, or Ge- 
ometric Optics approximations) in terms of \/C, where A 
and C are the wavelength and typical plasma length scale, 
respectively. In this approximation it is assumed that the 
electric and magnetic fields have a slowly varying amplitude 
with a rapidly varying phase, i.e. E^, _B M oc exp (iS) where 
S is the action, and V^S 1 = fc M defines the wave four- vector. 
Then, to first order in X/C, Maxwell's equations are 

fc M - E"u" + e"^ u a Bp) = W M E" , (7) 

k„ (Vb" - BV + e v " ais u a Ep)=0. (8) 



At this point it is useful to point out a number of proper- 
ties of E^ and _B M that follow directly from their definitions 
and Maxwell's equations. 

(i) u^E^ — u^B 11 = 0, which follows directly from the 
definitions of _E M and B 11 and the antisymmetry of F^ v and 
"F^" . 

(ii) = 0, which follows from equation JSJ and the 
definition of B M . 

(m) = 0, which follows from dE^B* = E li k v *F' iv = 

0, where u) = — fc M u M (chosen so that u) is positive) is the 
frequency in the LFCR frame and is assumed to be nonzero. 

(iv) ujB^B^ — — e 11 "™ 13 B^kvUcEp, which also follows from 
equation ©, luB^B^ + e^" 13 B^u a Ep = B^k v *F llu = 0. 

Properties (i)-(iv) define B M in terms of _E M , and w M : 

B" = --e» vaf) k v u a Ep. (9) 

LU 

Substituting equation (|5J into equations (J^J and 
gives 

LU 

*F^ = k a Ep . (If) 

LU 

Inserting these back into Maxwell's equations and combining 
yields 

Q^E" = , (12) 

where 

ft% = (k a k a 6Z - -Amwa^) (13) 

defines the dispersion tensor. 

Note that this is extremely general, all of the local 
physics is contained in the conductivity tensor. The expres- 
sions for the electromagnetic field tensor and its dual are for 
the radiation fields only. Hence, external fields appear only 
in the conductivity. 

2.2 Ray Equations 

Rays are well defined in the context of geometric optics. 
These are curves which are orthogonal at every point to the 
surfaces of constant phase (5). Given a relation in the form 
of equation 11211 it is possible to explicitly construct these 
rays . This has been done in detail for Euclidean spaces (see 
e.g. IWeinberdlT962Tl . The generalisation to a Riemannian 
space is s traightforward and will be done in analogy with 
IWeinberel <|l962h . 

Consider the general case of an equation governing the 
dynamics of a field, $ , in space time in terms of a linear 
operator, M, 

M(V M ,x M )* = 0. (14) 

Expanding in a two length scale approximation, as in §2.1, 
gives to lowest order 

M(fc M ,a; M )* = 0. (15) 

This implies that det M (fc M , x^) — along the rays of the 
wave field. This provides a dispersion relation, D (fc M , x^), 
a scalar function of the wave four-vector and position that 
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vanishes along the ray. If the eigenvalues of M are nonde- 
generate, then this also uniquely defines the polarisation of 

The ray can now be explicitly constructed by employ- 
ing the least action principle. The action can be explicitly 
constructed from the wave four-vector and the position by 



5(ri,T2) = 



fc, 



dx" 
' dr 



dr , 



(16) 



where r is an affine parameter along the ray. Let V be the 
hypersurface passing through the point x" (ri). By defini- 
tion, k u (n) is perpendicular to V. By varying S (ri, ts) with 
respect to fe M and restricting i M (n) to lie on V, it is pos- 
sible to derive equations which define the ray, 



5S 



T2 



^Sx 

dx" 

dfe M dx" 
dx" dr 



, dx" 

~~d7 



dkfj, dx" 
dr dx v 



dx" 
~dr~ 



dr 



Sx"dr + kn 8x"\ 



k^Sx"] 



(17) 



0. 



Because x" (ri) is restricted to lie upon V 
Because at T2 it is necessary for fc M (T2) = V U S — > 5S 1 = 
k/j, (T2) (Jx' 1 (T2). These imply that the integral must vanish 
for arbitrary variations. This will be generally true if 



dx" _ 
~d~r~ ~ 

and hence, 

dk u 



dP 

dk^ 

dx" 



and 



dk^ 
~dT 



(18) 



dx" dr 



dk^ dx" 
dr dx v 
dP \ dk u 



dk u 



dx' 



+ 



dP 

dxv- 



dx" 
dx" 



= 0, 



_ dP dr 
dr dx" 

where the final equality follows from the fact that P is 
constant along the path (namely P(k l _ ll x") — 0). There- 
fore, equations (11811 can be used to construct a ray given 
initial conditions and a dispersion relation. These are co- 
variant analogues of Hamilton's equations. Note that the 
affine parameterisation depends upon the particular form 
of the dispersion relation. For example, from D' (k^jx") = 
f (fc M , x") P (k u , x") it is possible to construct the rays as- 
sociated with P = 0, with the affine parameters related by 
dr' — dr/f: i.e. 

dx" 
~d~? 



f 



dP' 
dk u 

dP 

dk u 



+ P 



df_ 

dk u 



~ 1 dr ' 



(19) 



and similarly for fc M . Hence, any convenient affine param- 
eterisation can be selected by employing the appropriate 
function /. 

While this derivation is done in some generality, in this 
paper M = Q," v and * = E" . 



be obtained via kinetic theory. Three assumptions are made 
in the derivations below; (i) the equations of motion of the 
electrons are well approximated by the lowest order pertur- 
bations, (ii) the motions of the electrons are non-relativistic, 
and (Hi) the electrons execute motions over a small enough 
region of space that all other forces may be considered con- 
stant. Assumptions (i) and (ii) are often employed in stan- 
dard plasma physics. Assumption (Hi) will generally be true 
as long as the geometric optics approximation holds. 



2.3.1 Isotropic Cold Electron Plasma 

This is considered as an example and a limit of the case 
where a con stant external magnetic field is applied (cf. 
lDendvlll990l) . 

It is useful to introduce an order parameter (e) to lin- 
earise the force equations. All field quantities are clearly 
of first order. In addition, the change in the velocity of 
the charged particles is of first order (bu" = u" — u" oc 
eexp(iS')). Then, the electromagnetic force upon a single 
electron is given by 



ttm _ f" v P „ 

J EM — ' HU/lt 



(20) 



— u ml* ma — v 1 ufafl — j-> 

-eu eE u v — eeE u u v + ee u a eBpu v . 



In the first and third terms only the deviation from u" 
contributes, thus they are of order e 2 . In the second term 
u"u^ — —1 + 0(e) hence there is a first order contribution, 
and fjj, = eE". The force is related to u" to first order 
in e by J~em = —ium 5a" . The current is related to bu" by 
J" = enSu" . Therefore, the conductivity tensor is given by 



Up 



■ -8% , 
47riu " ' 



(21) 



where up = ^/4ire 2 n/m is the plasma frequency. 



2.3.2 Magnetoactive Cold Electron Plasma 

In the presence of an externally generated magnetic field, 
B" , (defined in the LFCR frame in the same way as B"), 
the electromagnetic force upon a single electron is 

^ M =F"»eu„ (22) 
= eu" eE" u v — eeE"u L ' u v + ee" VOi ^ u a (eBp + Bp) u v . 

In contrast to equation (I2(H . there is a first order contribu- 
tion from the third term in this case. Hence, to first order 
J"em = eE" + ee'""*' 3 u a Bp u„. It is useful to decompose bu" 
and E" into temporal, and spatial components along and 
orthogonal to B": 



5u^(5u v u")u" , hf=SlB" 
5uj_ = 8a" — bu" — Suu , 



(23) 



2.3 Ohm's Law for Cold Plasmas 



At this point it is necessary to determine an explicit form 
for the conductivity tensor a" v . For cold plasmas this can 



EC 



B V E V 
B a B a 



E" 



E" -EC 



(24) 
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With these new definitions it is simple to show that the force 
equation separates into 



-iujduf — , 

11 m 1 

-itvdu^ = — E'f H e w a y; a J3„ ftu^ v i 

m m 



(25) 



Clearly Jy 1 = — (up/^Tviuj^ E^. The perpendicular compo- 
nent may be determined by taking a second proper time 
derivative whence, to lowest order, 

— w 2 &tf = — iuo — E^ 
m 

+ —e' ll ' al3 u a Bp \—E± v + —E viaE u 7 B''5u 1 , 
m V m m 

= -iu;±E f [+(-y s flva0 u a B l 3E ±v (26) 
m \mJ 



— B V B V &tf . 

Defining oj s = (e/m) 2 B M S M and solving for = en&tj 
gives 



<5uf = 



2 



wsT + ^-er^uME^. (27) 

47T (dJg — (iHj V m 



After substituting in the expressions for i?^ and E^ the total 
current is given by 



2 



(tjfj — Ul 2 ) 



2 uk . 2 



(28) 



m 



As a result, the conductivity tensor can be identified as 



2 

Ulp 



>nv — — A — : — 7— 2 5\ I ~ w 9v- v + 

Amu; (ujg — 



ByB^_ 

B a B a 



ico—e^islf B p ) . (29) 
m 



In a flat space, the spatial components of this can be com- 
pared to the stan dard result (see e.g. iBovd fc Sanderson! 
ll969l:lDendvlll99d) . 



2.4 Ohm's Law for Warm Plasmas 

For AGN and X-ray binaries, accreting plasma near the cen- 
tral compact object will in general be hot. Even in low lumi- 
nosity AGN, acc reting electrons can have Vs on the order o f 
10- 10 3 (see e.g. lMeiia fc Falckel200ll : lNaravan et al.ll99H . 
In these environments assumption (m) in H2.3I that the mo- 
tions of the electrons are non-relativistic, is no longer valid. 

For warm plasmas, ones in which the thermal ve- 
locities of the electrons are significant compared to the 
phase velocities of the modes, it is possible to determine 
the conductivities using the Vlasov equation just as in 
flat space (see e.g. lDendvlll99d: IBovd fc Sanderson! Il969t 
iMontgomerv fc Tidmanlll964l) : 



df_ 



df_ 



0. 



(30) 



pi' \ r / xl L 

where and / are the momentum and distribution function 



of the electrons, respectively. The average plasma velocity, 
u 11 , must now be averaged over temperature in addition to 
the induced oscillations. Note that unlike the analyses of 
warm plasmas in flat space, this must now be done in a 
manifestly covariant way. At this point it is necessary to de- 
termine the form of the force, ?e M , under which the system 
is evolving. 



2.4-1 Isotropic Warm Electron Plasma 

In this case J-£ m = F^ v eu v . Hence expanding the distribu- 
tion function in terms of the order parameter introduced in 
H2.3.1l to first order, f = fo + e/i + 0(e ), and inserting into 
equation (130ft gives 



(31) 



Considering the lowest order in the two length scale expan- 
sion of H2.ll this may now be solved for fj : 



_ ieu v f dfo 
~ u a k a \dp>* 



(32) 



which is the covariant analogue of t he expressions found in 
the kinetic theory literature (see e.g. lDendvlll990t) . 

Assuming that the plasma was originally charge neu- 
tral the current density is related to the perturbation in the 
distribution function, /i , by 

J M = e / d 4 p /iu M . 



Then, using equation 1101 this may be written in terms of 
£ M as 



r = - — k a E v 



4 U 

a p 



dfo 
dp v 



dfo_ 

dp a 



(33) 



From this it is clear that the conductivity tensor is 



v u = k / dp—— 

urni J p p Kp 



Pa 



dfo 
dp" 



Vv 



dfo 
dp a 



(34) 



In order to make a connection with the expression derived 
in the previous section it is convenient to integrate this by 
parts, 



9 in, 



p a k a 



+ 



fo , (35) 



where the boundary terms vanish by virtue of the conver- 
gence of J d*p fo- For the cold plasma, fo = nS\p^ — mu^), 
thus, 

cop ( kuV-v ~\~ kviiu k a k Un l Uv\ /or.\ 
a ^=-^( 9 ^ + u + — — • (36) 



This differs from the result in H2.3.1l in two respects: terms 
proportional to and the term proportional to Because 
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the conductivity enters Maxwell's equations only through a 
contraction with the electric four-vector, the former are su- 
perfluous. The latter represents the sonic mode which ap- 
pears in the kinetic calculation of the conductivity only in 
the form of an infinite wavelength mode. For the two trans- 
verse electromagnetic modes (E^k^ — 0) this does agree. 



2.4-2 Magnetoactive Warm Electron Plasma 

In the presence of an external magnetic field has a 

zeroth order contribution: 



■J F. 



eF^ u v + eF£" it„ . 



(37) 



where, in terms of the external magnetic field (again defined 
in the LFCR frame), F£ = E-^^uJSp, (cf. equation Q). 
Expanding the Vlasov equation in the perturbation param- 
eter e to first order and in the two length scale expansion 
(j EJ now gives, 



m \ op^ 



=-±F^m 

m \ op^ 



(38) 

At this point it is useful to introduce a function n defined 
implicitly by 



± = ^F» v p v ( — 
dr] m Ex " V dp' 1 , xfi 



(39) 



fcf. lLifshitz fc Pitaevskill98ltlKrall fc Trivelpiecei 19731) . In 

terms of 77, the electron momenta are determined by the 
equation 



dp* 1 e. y-^ui/ c uvctft — 

—r- = —F£ x p v = — e M v UaBpPv ■ 
dr] m m 



(40) 



As in the cold case, this may be reduced to a two dimensional 
problem by an appropriate decomposition of the momentum: 



LI Lt LL Lt 

P± = P ~ Pt ~ P\\ i 



(41) 



In terms of these, the system of equations for p^ reduce to 

dtf_ 
di] 

dpi 
dr] 



= 0, 
= 0, 



(42) 



u r| l " LLUap YD 

dr] m 



This last equation is simply that governing cyclotron mo- 
tion. Using the fact that d/dr/ commutes with the metric 
(this is because the metric depends only upon x^ and not 
p M ) it may be rewritten as a pair of uncoupled, second order 
ordinary differential equations: 



d Px , 2 M 



0. 



This has solutions 



= p£ cos(u>bV + <j>o) + Py sm(cu B r] + <j>o) 



(43) 



(44) 



where p£ and pt^ are a pair of bases which span the space 
perpendicular to u M and £> M , and <f>o is a phase factor. By 



inserting this solution into equation 1401 and matching up 
trigonometric terms, p£ can be found in terms of p£, 



Py = — £ UaBaPxp. 

0JB 



(45) 



It is possible to now solve for 77 in terms of p M , p%, and <j>o: 



V =± arctan e£ ^ M y agd . (46) 

Ub [ y muaPxPt J 

Inserting p M (»y) into /1 and fo transform equation 138> 
into a first order differential equation for f\. This has solu- 
tion 



h 



/' 



1 J p/V 



dr] E" 



where 



dfo 



dfo 



H = exp 



%k u 



-dr) 



(47) 



(48) 



(49) 



The integral for /1 may be rewritten in terms of p M by using 
equations 1421 and 1431 1 . 



J (pt + pfj dr] = [pi + p^j 77 , 
d 2 pl 



P± d V = —2 



drf ujg dr] 



- VLVOiQ k 1 

— £ U a Df3p±v. 



(50) 



(51) 



Thus, 

]i = exp< 1 



B v k v 



B a B c 



8m )V 



2 fe W 



(52) 



With equation 14411 this may be treated as a function of 77, 
while with equation 1461 this may be treated as a function 
of p". 

As in the previous case, the current four-vector is then 
found by integrating over the momentum portion of the 
phase space. This gives the conductivity tensor to be 



e 

777 



d 4 pp M 



(53) 



where it has been emphasised that the interior integral is to 
be treated as a function of the momenta. 



2.4-3 Conductivity in Quasi- Longitudinal Approximation 

In general, the integrals over 77 in equation I53H can be 
evaluated in terms of sums of Bessel functions in an anal- 
ogous fashion to that typically done for the non-relativistic 
case (see e.g. iKrall fc Trivelpiecd H97I . Nonetheless, this 
can be significantly simplified by considering the case where 
(i) fo is a function of V 2 = p m Pm aQ d £ = p M u M only 
(typically fo can be written in the form f(e)S(V 2 — m 2 ) 
where the delta function is required to place the distribution 
on the mass-shell), («) e ,lva ' 13 TicBpk^ — (i.e. the quasi- 
longitudinal approximation), (lit) u>b <C U), and (iv) fo is 
such that p^u^/m -1<1 (i.e. cool, not hot). 
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Assumption (i) simplifies f) t 



^ @J0 j / U V V Ll\ 

7T- % (U P ~ U P^) 

mu> de 



(54) 



Note that because e is independent of rj, the terms involving 
fo can now be brought out of the innermost integral in equa- 
tion 1531 . Assumption (n) gives that = and hence, 



These may be integrated by parts to produce 

JM" _ 1 IK i: I ./<■' ' 



1 2 



d*P fo [ s M " 

,4 fo 

dp — 



mzu 
2p"k" + p v k tl 



+ ft K a 



(64) 



fj, = e 



(55) 



where = k^p^/m. Therefore, the two integrals that must 
be done are 



p|7 e dr\ — p. 



(56) 



and 



/ 



u a Bp ) p±„ 



/'• 



(57) 



Therefore, in the quasi-longitudinal approximation, 
e dfo 1 



fi = 



iwm de zu 2 



2 2 B^Bv 



3U B 



p v E" 



(58) 



where the definitions of pff, p^_, and _E M were used. In the 
quasi-longitudinal approximation, is orthogonal to the 
external magnetic field, B^. As a result, the there are only 
two integrals that must be done in order to find the conduc- 
tivity tensor: 

iuj f 4 izu p v dfo 
I « V — — V 

m 



, — TP P a 



^2 = / dp — 



2 PP 



9e 



In terms of these, the conductivity is 

„2 



e 

itom 



(59) 



(60) 



(61) 

u; ae 

Noting that the will only be contracted on the second 
index with terms orthogonal to « M (for If this is the electric 
field), the I^ v are given by, 

p»k a \ dfo 



From equation 1541 it follows that 

dfo _ Pft k a / df dfo 



If = -i / d 4 p 



-yr w 



m J dp 



*2 = / d P—> TP (™9 ~P k -7T— ■ (62) 

m J zu 2 — u 2 B dp a 

Because there is already a term linear inwg in equation 
I6UII . to lowest order in assumption (in) uj% may be neglected 
in the Thus, 

4 n f ua P"k a \ dfo 



I£ = / dpp* 1 g 



l T = -* / d P — 



mvj J dp a 

P" ka \ d J±. ( 63 ) 



mzu J dp' 



Note that in this case, is simply the integral that had to 
be done for the warm isotropic plasma (cf. equation J35J). 

Assumption (iv) enters by expanding zu about u>. De- 
fine p 2 = e 2 — m 2 , i.e. p is the magnitude of the spatial 
components of the momentum in the LFCR frame. Then, to 
second order in p, 



(-cuf 



1- j 



P^B^B, 



muB a B a 

,/(./- 1) ( p-B.I. : B. 

mujB a B a 



+ 3 



2m? 



(65) 



Thus, 

if" 
1 i 



J 2 



d 4 p fo 



uii . P P t R-o 

9 +—9. — -r 



fa 



2m : 



+ 



p»B^k v B v 
mujB a B a 

p^p" k a k, 



+ 2 



(66) 



where terms odd in p M and terms oc k v have been dropped. 
The former is due to the fact that fo has been chosen to 
be an isotropic function of the spatial components of the 
momentum in the LFCR frame and hence any odd terms 
will vanish upon integration. The latter is allowed because, 
as stated earlier, these will only have significance when con- 
tracted with terms orthogonal to k^ (for 7^" this is results 
from the quasi-longitudinal approximation in which fc M can 
be written in terms of u M and £> M only). From symmetry it 
is clear that 



4 t (P^m) 1 2 , t > 

dpfo ^B7 = r m </o> = 



where 



</o> 5 



j4 c 2 
dp fop ■ 



(67) 



(68) 



In addition, the off-diagonal components of the integrals over 
p^p v will vanish due to the symmetry of fo- Because adding 
terms oc tt will not alter the physical solutions, it is possi- 
ble to replace J dSpp^p" fo with |nm 2 (/o) 2 '• Lastly, note 
that 



(k v B v ) A 



B a B a 

Therefore, the are given by 
If ~ nT^ v and 



uj l + k a k a 



1 2 



(69) 



(70) 



where 



1 k a k 

T lS l + I^(/ ) 2 

l2 = l-\ </0> 2 + ^ 

6 w J 



;/o) 2 



(71) 
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Because the terms multiplying I2 in the conductivity are 
already of first order (the order of ujb is necessarily equal to 
or smaller than that of p for the approximations thus far to 
hold) , to second order in small quantities in the conductivity, 
T2 — 1. As a result, with the lowest order finite temperature 
corrections the conductivity is given by 



UJp 

4:TViUJ 



(72) 



For the cold plasma T\ = 1 and this does reduce to the 
appropriate expansion of the conductivity derived in H2.3.2I 



2.5 Dispersion Relations 

Given the conductivities derived in H2.3I & H2.4I it is now 
possible to obtain the associated dispersion relations. It is 
instructive to compare these to the dispersion relation for 
massive particles (de Broglie waves): 



Dfkift , x ) — k kfj, -\- 



(73) 



That this does produce the time- like geodesies when inserted 
into the ray equations is demonstrated in appendix 1X1 



2.5.1 Isotropic Electron Plasma 



The conductivity tensor obtained in ^2.3. II for the isotropic 
cold electron plasma yields the dispersion tensor 

= (k a k a + u%) - fc'X . (74) 

For the transverse modes, this gives the dispersion relation 

£>(*„, z M ) = *"* M + o£, (75) 

(cf. iKulsrud fc Loeblll992l) . For constant density plasmas 
this is nothing more than the massive particle equation, 
cf. equation l|73^. For plasmas with spatially varying den- 
sities this leads to a variable effective "mass" . Hence in gen- 
eral, photons in plasmas will not follow geodesies. This is a 
representation of the refractive nature of the plasma. 



2.5.2 Quasi-Longitudinal Approximation for the Cold 
Electron Plasma 

When magnetic fields are present it is necessary to utilise 
the conductivity tensor obtained in H2.3.2I In the quasi- 
longitudinal approximation the wave four-vector is parallel 
to the external magnetic field. In this approximation, the 
modes are transverse. This follows from the fact that in the 
LFCR frame this is true and that since this is a local prop- 
erty expressible in covariant form, it must also be true in 
an arbitrary frame. This can be explicitly verified by com- 
parison with the results of H2. 5. 51 where the general case is 
considered. 

Under these conditions the dispersion tensor takes the 

form 



where a, 7, and M M „ are defined by 



(76) 



a = k^ku — 5uj , 7 = Suj ( — ) , (77) 



2 

UJp 



2 j ±v*fj.v 



Taking the determinant of ^V L V yields 
det fit = a - a 2 7 2 B M B M 



= a 2 (a — Sujlub) {p. + Sloujb) = . 



(78) 



The two modes corresponding to a — are the sonic mode 
and the unphysical mode proportional to u fl which is elimi- 
nated by the condition that u^E^ = 0. The other two modes 
have dispersion relations 



D (k/j,,x^) = a ± Sujlob 



k^kn -)- 



UJUlp 



(79) 



LO ± UlB 

As with equation 17511 . this dispersion relation also has a 
term that could be identified with the mass in equation 
1731 1. In contrast with equation 1751 . now that "mass" de- 
pends upon the polarisation eigenmode. As a result, differ- 
ent eigenmodes will propagate differently. Again this is an 
expression of the dispersive nature of a magnetised plasma. 

In addition to dispersion, a noticeable departure from 
its non-relativistic analogue is the presence of fc M in the defi- 
nition of ui. This is not surprising since it is the most general 
Lorentz covariant extension of the quasi-longitudinal disper- 
sion relation. Of interest is the fact that the dispersion re- 
lation is now cubic in the magnitude of k, k. Because two 
roots clearly exist in the low density limit, a third root must 
also exist. This results in a new branch in the dispersion 
relation. This will be explored in more detail in CHI 



2.5.3 Quasi-Longitudinal Approximation for the Warm 
Electron Plasma 

For the conductivity derived in H2.4.3I this is identical to the 
previous section, where a and 5, are replaced by k^k^+Tiojp 
and —Lup/uj 2 . Then, 

D (ka , £ M ) = ct ± ujp 

UJ 

1 a , . t- 2 1 2 Wfl 

= k ku + l\u)p ± ujp 

UJ 



= 1 + 



1 Up 

37 



</o) 2 J fe^ + wjiw 2 ,^. (80) 



For a thermal electron distribution, (/o) 2 — 3fcT/m and 
hence 



1 ujp u4 

\/°>2 = ~2 where 



2 _ kT 2 

UJip — OJp . 

m 



(81) 



Note that lot is related to the Debye frequency, lod, by ujt = 
ujp/uJD- Thus, including the lowest order finite temperature 
corrections, the dispersion relation in the quasi-longitudinal 
approximation is 



£>(fc M ,0 



1 + 4)^ + ^*4-. (82) 
uj J UJ 



2.5.4 General Magnetoactive Cold Pair Plasma 

The conductivity for the pair plasma may be obtained by 
adding the conductivities for the electrons and the positrons, 



0~n.11 = 0",, 



+ 0~ \iv 
2 

UJp 



A-Kiu) (ujp — uj 2 ) 



-v ^+^ BWK ) , (83) 
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where now the plasma frequency is defined in terms of the 
sum of the number densities of the electrons and positrons. 
The resulting dispersion tensor is 



Q,^ = ag^ v — k^k v + (3B^B„ . 



(84) 



where a, 7, 5, and M MI/ are defined as in equation 1771 . 
and f3 = 5{e/m) 2 . In addition to the requirement that 
Q.^ V E V — 0, must be orthogonal to w M . As a result, it 
is necessary to alter fi M „ in such a way that it explicitly sep- 
arates the eigenmodes orthogonal to u M from the unphysical 
mode. This can be trivially accomplished by adding a term 
—cokf/tit, to the dispersion tensor. Note that this does not 
change the dispersion equation for the physical modes be- 
cause E^u., — 0. Thus, consider 



MI/ = ag^ v — fc M (k v — uu v ) + pB^Bv , 



(85) 



instead of the dispersion tensor given in equation 184^ . For 
this dispersion tensor, the unphysical mode is trivially found 
to be , with dispersion relation D — a. As in ^2.5.21 the 
dispersion relations can be found by taking the determinant 
of the dispersion tensor: 



detf^ 



-(1 + <5)wV 



a + 5uj b 
5 



1 + 6 



eB»k L 



(86) 



where the definition of a was used. Therefore, the dispersion 
relations for the two electromagnetic modes are 



D 1 (k„,x tl ) = fc M fc M - 



Up 



(87) 



2 

LOp 



Ulp + LOg 



eB^k L 



It is straightforward to show that Di and D2 correspond 
to the extraordinary and ordinary modes, respectively, by 
considering the transverse limit (B^k^ — 0). 

2.5.5 General Magnetoactive Cold Electron Plasma 

For the general case, no approximations, except those used 
to derive equations 1121 and 1291 , are made. In this case, in- 
serting the conductivity tensor obtained in H2.3.2l into equa- 
tion M'6l gives 

Qiiv = k k a g, LV kiikv - ~ — i io gnu 

(«B ~ "> ) V 

Collecting the coefficients of like tensors gives 

Q M „ = as", - fc M fc„ + (5B"B V - i'yMl , (89) 



where a, (3, 7, 8, and M^u are defined as in >!2.5.2l and ^2.5.41 
As in the previous section, it is useful to add a term propor- 
tional to u v E v to the dispersion equation. Hence consider 



■■ - fc M (k v - uu u ) + fiB^Bu - i-yM* . 



(90) 



Proceeding as in the previous sections, the scalar dispersion 
relations corresponding to the different eigenmodes can be 



found by considering the determinant of the dispersion ten- 
sor: 

detfi M „ = aja 3 + [pB^Bp - (fc"fc M + uj 2 )] a 2 
— 5lo 2 b (fc M fc M + uo 2 ) — 5 ^— £> M fc M j + 8 2 ijj 2 lo% 



8uj% - —B' l k 



(91) 



Inserting the definition of a reduces the terms in the braces 
to a quadratic in k^k^, which may be solved to produce the 
desired dispersion relation: 



k^ku — 8u} 



2(1 + 5) 



eB^k u 



(l + 2<5)wl 



eB^kaV 



( eB^k. 



■ 2 (2a; 2 - u% - wf 



(92) 

This is a covariant exte nsion of the Appleton-H artree dis- 
persion relation (see e.g. lBovd fc Sandersonll 969'l. As in the 
previous two sections, this continues to bear a resemblance 
to the dispersion relation for massive particles. Again the ef- 
fective "mass" depends upon position and the polarisation 
eigenmode. Additionally, it now depends upon the direction 
of propagation relative to the external magnetic field as well. 



3 EXAMPLE APPLICATIONS 

In 35] the general theory of a covariant magnetoionic the- 
ory was presented for electron-ion (in the Appelton-Hartree 
limit) and pair plasmas. While astrophysical plasmas will in 
general be warm, the cold electron plasma does provide a 
instructive setting in which to highlight some of the similar- 
ities and differences that a fully general relativistic magne- 
toionic theory has compared to general relativity or plasma 
effects alone. 



3.1 Bulk Plasma Flows 

A number of effects will appear in special relativistic plasma 
flows. The covariant formulation of magnetoionic theory can 
have implications for the structure of the dispersion relation. 
As briefly mentioned in H2.5.2I the equation for the mag- 
nitude of the spatial part of the wave vector is now cubic. 
This is essentially due to Doppler shifting. Thus these effects 
should appear in relativistic bulk plasma flows as well as in 
regions of strong frame dragging (e.g. near the ergosphere of 
a Kerr hole). 

For a relativistic bulk flow (in the x direction) 

k t - vk x cos 8 
U= Vl-v 2 ' (93) 
where 9 is the angle between the wave vector and the motion 
and v is the velocity of the motion. Clearly the coupling be- 
tween the previously mentioned third branch depends upon 
both v and 6, being strongest when — 0. Shown in figure 
0are the quasi-longitudinal dispersion relations for a rela- 
tivistic bulk flow for a number of velocities and magnetic 



© 0000 RAS, MNRAS 000, 000-000 



Covariant Magnetoionic Theory I: Ray Propagation 9 




Figure 1. The dispersion diagram at a number of magnetic field 
strengths and velocities for a relativistic bulk plasma flow. The 
frequency scale is set by = 1. The ordinary (extraordinary) 
eigenmode is shown by the thick (thin) line. Note that the dis- 
persion diagrams are asymmetric due to the plasma motion. 



field strengths and 6 = 0. The frequencies are measured in 
units of the plasma frequency, making this otherwise scale 
invariant. Note that a whistler-like branch appears for the 
ordinary mode which is not present in the non-relativistic 
theories. Similar to the whistler branch of the extraordinary 
mode, it is asymmetric due to the bulk motion. In the limit 
of vanishing plasma density this branch does not transform 
into a vacuum branch, in much the same manner as portions 
of the whistler. Therefore, this mode cannot escape from the 
plasma, necessarily reflecting at the surfaces of the plasma 
distribution. This may have implications for the pressure 
balance in thick disks with large velocity shears and jets, 
even at frequencies where these are optically thin. 

In bulk plasma flows the new branch appears because 
the velocity mixes the spatial and temporal components of 
fc M . In a Kerr spacetime, frame dragging is responsible for 
mixing these components. In this case 



" fee 



(94) 



This is similar to equation I93H with the role of the velocity 
being taken by g^ /g U - Hence, the overall effect is qualita- 
tively the same; a new branch similar to the whistler appears 
for the ordinary mode. 



3.2 Isotropic Plasmas and Particle Dynamics 

In both special and general relativistic settings, the prop- 
agation of photons through an isotropic (field free) plasma 
can be represented in a manner ana logous to that of parti- 
cle dynamics in a potential (see e.g. iThompson et al. 1994. 
for the non-relativistic case). Following the manipulations 



in appendix [X] it is straight forward to show that for the 
dispersion relation given in H2.5.1I D — k^k^ + up, that 



-V m 4l^ 



(95) 



where u M = dx^/dr, i.e. 4wf> acts as a potential in which 
the the photons propagate (the factor of 4 is due to the 
particular affine parameter chosen, namely that associated 
with the choice of the dispersion relation given above). 

For plasmas in which magnetoionic effects are not signif- 
icant to the photon propagation (magnetoionic effects may 
still be important for emission and the propagation of po- 
larisation) this allows a somewhat more simplified analysis. 
If enough symmetries are present, then the rays may be de- 
termined via direct integration. For example, consider a sta- 
tionary, spherically symmetric plasma distribution around a 
Schwarzschild black hole. In this case equation Ij95^ shows 
that vt and Vs are conserved, associated with the time and 
azimuthal Killing vector fields, respectively. Therefore, with 
the dispersion relation, 



dt 



= v — g v t = 



1 - 



2M 



I't 



dr r 



L , , 



(96) 



Which may be directly integrated to give the ray as a func- 
tion of the affine parameter r in precisely the same fash- 
ion as is typically done to find the particle orbits of the 
Schwarzschild metric. 



3.3 Photon Capture Cross Sections 

In the vicinity of a black hole, polarisation can arise even in 
the case of a gray emissivity. This occurs when one mode is 
preferentially captured by the black hole due to dispersive 
plasma effects. Even without a method for performing the 
radiative transfer, this can be estimated by considering the 
photon capture cross section of Schwarzschild black hole. 
It is necessary to provide a plasma geometry - the plasma 
density, velocity, and magnetic field - as functions of po- 
sition. Here, the density is given by the self-similar Bondi 
solution, tup oc r -3//4 . The magnetic field is chosen to be a 
fixed fraction of the equipartion value, ljb oc r~ 5 ^ 4 . Finally, 
the velocity is chosen such that the plasma has zero angular 
momentum, i.e. Tit = 9 U an d u r = ug =u,p = 0. While 

this doesn't correspond to a realistic accretion flow, it does 
provide insight into the type of effects dispersion can have. In 
order to further simplify the problem the quasi-longitudinal 
approximation was used. Typically this is a good approx- 
imation, only failing when the angle between k^ and S M 
is within ~ upuip/^ of ir/2. This dispersive polarisation 
mechanism produces primarily circular polarisation for the 
same reason. 

Shown in figure |2] are these cross sections for a num- 
ber of different plasma densities (through wp) and magnetic 
field strengths (through ujb)- These are both scaled by the 
observed frequency at infinity, and hence are not tied to 
any particular frequency scale. The capture cross section of 
the extraordinary mode decreases more rapidly than that 
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Figure 2. Photon capture cross sections in units of the vacuum 
capture cross section, <r 7 = 27ttM 2 , for the quasi-longitudinal ap- 
proximation as a function of plasma density (u!p/u> ohB is the value 
of the plasma frequency at r = 3M) at a number of magnetic field 
strengths. The solid, dotted, short dashed, long dashed, and dash- 
dotted lines correspond to u>p/u> obB = 0, 0.7, 1.4, 2.1, and 2.8, re- 
spectively, at r = 3M. The inset shows the circular polarisation 
fraction, m c , in terms of the effective emission area, A, for the 
same set of magnetic field strengths. 

of the ordinary mode, with increasing density. The dispar- 
ity between the two capture cross sections increases with 
increasing magnetic field strength. 

This can be a very efficient manner of creating polar- 
isation over the inner portions of the accretion flow. How- 
ever, far from the hole (outside the inner 5 — 10M) this 
becomes a small effect. As a result, the fraction of polarisa- 
tion produced depends upon the magnitude of the diluting 
emission from regions of the accretion flow distant from the 
hole. Nonetheless, it is possible to parameterise the unknown 
emission in terms of an effective emitting area (the details 
of which still depend upon the details of the accretion flow). 
Shown in the inset of figure is the circular polarisation 
fraction scaled by the effective emission area in units of the 
vacuum photon capture cross section. 

3.4 Tracing Rays 

With general dispersion relation for cold magnetoactive plas- 
mas, equation 1921 . and the ray equations, equations list , it 
is straightforward to explicitly construct rays. The plasma 
geometry outlined in the previous section will be used here 
as well, with the scales set by uip(r — 3M) = aj obs and 
ub{t = 3Af) = 2uj abB , where o; oba is the frequency observed 
at infinity. In figure [3] rays are propagated in the vicinity of 
a Schwarzschild black hole. For comparison, in figure 0] rays 
are propagated near a maximally rotating Kerr black hole. 
The null geodesies are shown by the dotted lines for refer- 
ence. In both figures the extraordinary mode (solid lines) is 




-20 -10 10 20 



y(M) 

Figure 3. The paths of the ordinary and extraordinary polarisa- 
tion eigenmodes in the vicinity of a Schwarzschild black hole are 
shown by the dashed and solid lines, respectively, for a number of 
impact parameters. The dotted lines show the null geodesies for 
comparison. The x axis lies along the ray paths at infinity, and 
the y axis is orthogonal to both, the x axis and the slice of im- 
pact parameters considered. The plasma density is ex r -3 / 2 and 
u)p(r = 3M) = u) obB . The magnetic field has a split monopole 
geometry with its strength oc r~ 5 / 4 and ug(r = 3M) = 2u; b s . 
The horizon is shown by the filled region in the centre. 

refracted the most, and the ordinary mode (dashed lines) 
is refracted more than the null geodesies. This is precisely 
what is expected on the basis of the capture cross sections 
presented in H3.3I In addition to dispersive plasma effects, 
comparison with the null geodesies demonstrates that gen- 
eral relativistic effects are also significant. 

3.5 Intensity and Polarisation Maps 

The impact that dispersive plasma effects can have upon the 
spectrum of an accreting object can be illustrated by maps 
of the intensity. Here, in addition to the plasma geometry 
employed in the previous two sections, an optically thick 
Shakura-Sunyaev disk is introduced. The emission is solely 
from this disk and assumed to be thermal with 

/ / \ 3 / 10 

T(r) VR^jr) 

(see e.g. iFrank et al1ll992l) . The overall constant is depen- 
dent upon a number of disk parameters and hence is not of 
particular interest here. Nonetheless, it is chosen such that 
kT(r — oo) = v obB for convenience. The innermost radius of 
the disk, J? min , is chosen to be 3M. Doppler effects due to 
the rotation of the disk are ignored here. 

Shown in figure |^| are the intensity maps for when 
(a) plasma effects are neglected, (b) plasma effects are in- 
cluded, (c) only the left-handed circular polarisation (ordi- 
nary mode) is considered, (d) only the right-handed circu- 
lar polarisation (extraordinary mode) is considered. Because 
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Figure 4. The paths of the ordinary and extraordinary po- 
larisation eigenmodes in the vicinity of maximally rotating Kerr 
black hole are shown by the dashed and solid lines, respectively, 
for a number of impact parameters. The dotted lines show null 
geodesies for comparison. The plasma parameters are the same 
as those for figure [3] In addition to the horizon, the ergosphcre 
is shown by the partially shaded region. The rays originate from 
60° above the equatorial plane. The y axis is orthogonal to the 
the rotation axis of the black hole. 

the overall flux from the disk is dependent upon the details 
of the accretion flow, the intensities are normalised by the 
highest intensity in panel (b) . Comparing panels (a) and (b) 
demonstrates that including dispersive plasma effects makes 
a significant difference. This difference originates primarily 
from contribution by the extraordinary mode shown in panel 
(d). 

As implied by figure [5] the shadow the black hole casts 
upon the extraordinary mode is less than that cast upon 
the ordinary mode, which is in turn less than that upon the 
null geodesies. In addition to the differences in the overall 
intensities, there is a substantial difference between the con- 
tributions from the two polarisations as seen by comparing 
panels (c) and (d). 



4 CONCLUSIONS 

The covariant magnetoionic theory developed here is distinct 
in many respects from the non-relativistic theory. Firstly, 
it qualitatively changes the topology of the dispersion re- 
lations, adding an entirely new branch, as shown in H3.ll 
Secondly, it allows the inclusion of gravitational lensing ef- 
fects, vital for application to compact accreting objects. In 
addition, as shown in H3.3l and H3.4I dispersion due to plasma 
effects can have a significant impact upon the propagation 
of photons in a dense plasma environment near a black hole. 
As demonstrated in H3.5I this will lead to a modification of 
the spectrum. As a result, studies which neglect dispersive 



Magnetoionic Theory I: Ray Propagation 11 



i — i 




-10 10 -10 10 



(c) (d) 

Figure 5. Shown is the normalised intensity for an optically 
thick, Shakura-Sunyaev disc around a Schwarzschild black hole 
when (a) plasma effects are neglected, (b) plasma effects are in- 
cluded, (c) only the left-handed circular polarisation (ordinary 
mode) is included, (d) only the right-handed circular polarisation 
(extraordinary mode) is included. Note the different scales for to- 
tal intensities ((a) and (b)) and the polarised intensities ((c) and 
(d)). The disk is inclined 60° relative to the line of sight. £ is 
parallel to the equatorial plane. r\ is in the line of sight-azimuthal 
axis plane. The overall scale is set by the choice of observation 
frequency and the parameters of the disk and hence are not rele- 
vant here. The plasma geometry is the same as that for figures 151 
and El 



plasma effects may be inappropriate when the observation 
frequencies are near the plasma and/or cyclotron frequen- 
cies. 

On the other hand, because plasma effects have the ca- 
pability of altering the spectrum, it is possible for the un- 
derlying plasma to be observationally probed using polarised 
flux measurements. For example, if the horizon of the black 
hole in the Galactic centre can be imaged (cf. iFalcke et alJ 
2000), observations of the polarisation map could easily dis- 
tinguish the nondispersive from the dispersive case, placing 
limits upon the local magnetic field strength and plasma 
density. Integrated values for the polarisation could yield 
useful information about the environments of other accret- 
ing systems, such as X-ray binaries and pulsars. Because 
these effects can be expected to be confined to the decade in 
frequency surrounding the plasma frequency, they should be 
easily distinguishable from the effects of different accretion 
models. 
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APPENDIX A: GEODESIC MOTION IN THE 
DISPERSION FORMALISM 

Given the dispersion relation in equation 1731 . 

D(k ti ,x' J ') = k^k^ — m 2 , 

and the ray equations list . 

dx" _ f 9D\ dk^ _ f dD\ 

dr-{dkj xtl and dr --\dx») hii > 

it is possible to derive the geodesic equation. The partial 
derivatives on the right side of the ray equations are 



dD 

dk u 



2k" 



(Al) 



and 



dD\ ( dk a kpg a ^ \ 3g^ 
dx^) kii ~\ dx» ) kll ~ "dx" 

= k ga@,fj, • 



Combining the ray equations gives 

d 2 x" _ 2 dk^_ _ dk v g» v 
dr 2 dr dr 

-2k dXa99 " V +2a^ dk » 



(A2) 



-kk?k a g >11 ' g/3 v , a + 2g" v k a k p ' g a p,v 
-4fe a jfc /3 iff' M ' {g a v£ + gpv, a - g a p,v) 

dx a dx 13 
dr dr a/3 ' 



(A3) 



where the definition of the Christoffel symbol was used, 
Y^g = \g" v (g av ,p + gp„, a - g a p,v). Collecting terms on the 
left produces the well known geodesic equation: 

d 2 x" dx a dx r * u 

+ - — -r^ = o, 
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